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ABSTRACT 

We present a new derivation of the (D{a) angular distribution of the outgoing g-quark in 
the production process e~^e~ — 7, Z — > qq{g). In our calculation, we express the three- 
particle phase-space integration of the gluon-bremsstrahlung process in terms of a general 
set of analytic integral solutions. A consistent treatment of the QCD one-loop corrections 
to the axial-vector current deserves special attention. This is relevant in the derivation of 
the forward-backward asymmetry predicted by the standard model. Finally, we provide 
the full analytical solutions for the differential rates in closed form and conclude with 
numerical estimates for bottom and top quark production. 
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1 Introduction 



The recent discovery of the top quark has intensified interest for radiative corrections in the 
Standard ModeL In particular, processes involving top pair production are very attractive 
for closer investigation due to the absence of fragmentation effects. For bottom quark 
production, the measurement of the forward-backward asymmetry from Z-decays allows 
for a very precise determination of the electroweak mixing angle and the vector /axial- 
vector couplings to the fermions. The high accuracy of these experiments as carried out 
by the LEP and SLC collaborations requires theoretical predictions beyond Born 
level in the perturbative expansion of electromagnetic and strong couplings. 

In this regard, the calculation of differential cross sections for the electroweak pro- 
duction of massive quarks from e~^e~ annihilation is of crucial importance. Radiative 
corrections to this process are in general dominated by QCD effects, which can be calcu- 
lated perturbatively at energies sufficiently above the production threshold of the quark 
pair. 

A first analytical treatment of radiative corrections to the angular distributions in 
heavy-quark production was carried out in Ref. 0] at energies where the e~^e~ annihilation 
process is dominated by pure photon exchange. At higher energies, the 7 — Z interference 
has to be included This leads to a considerable complication in the three-body phase- 
space (gluon bremsstrahlung off a massive quark pair). Only recently f^, expressions for 
the bremsstrahlung process have been obtained by including cuts on the energy of the 
emmitted gluon. 

In this article, we present a systematic derivation of the QCD one-loop corrections 
to the differential production rate for the annihilation process e^e~ — 7, Z — > qq{g). No 
approximations are made and the whole three-body phase-space is considered. The full 
dependence on the quark masses is kept so that the final analytical expressions apply to 
the full energy spectrum where perturbative QCD is valid. 

Our approach relies on a general set of integral solutions completely describing the 
complicated phase-space with two massive fermions and one vector particle. The method 
involves special integration techniques and was devised in Ref. 0. It was used in various 
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calculations at 0(a) such as the total longitudinal or beam-alignment polarization for 
quarks produced in e+e^ annihilation [0,^3. In the present work, we have been lead to 
an additional different class of phase-space integrals due to the more intricate angular 
structure of the differential calculation. Furthermore, we have used two different 75- 
prescriptions in the derivation of the C-odd part of the differential cross section. There, 
the corresponding traces contain an odd number of the Dirac 75 and one faces immediately 
the problem on how to treat it consistently within dimensional regularization. 

The article is organized as follows. In Section 2, we set up the basic formalism for 
the calculation of the differential cross section da{e'^e~ '~f,Z ^ qq)/dcos6 (the angle 
6 is the usual scattering of the g-quark with respect to the electron beam axis) including 
Born approximation and QCD virtual corrections. Special attention is paid to the gluonic 
vertex corrections within different regularization schemes. In Section 3 follows a detailed 
presentation of the kinematics of the massive qqg phase-space which directly gives the 
angular distributions for the 0(a) tree-graph contributions. Compact expressions for the 
real-gluon corrections are derived. The C-odd partial rate is calculated in two independent 
approaches yielding identical results. In Section 4, we present the full analytical results 
for the differential rates and conclude with an explicit numerical analysis for bottom and 
top quark production. The article is supplemented by two Appendices A and B where 
details on the renormalization of the axial-vector current and the massive three-body 
phase-space can be found. 



2 Basic Framework and Virtual Corrections 

The differential cross section for the production process e~^e~ qq is a binomial in cos^, 
6 being the scattering angle of the tagged quark. It is common to introduce the structure 
functions cr^, a^, and ap as follows 

4^— = |(l + cos2^)cr^ + f sin^^d^ + f cos^a^, (1) 

so that au and cr^ are the contributions stemming from unpolarized and longitudinally 
polarized gauge bosons, respectively. The structure function ap relates to the difference 
of left- and right-chiral polarizations of the quarks and constitutes the C-odd component 



under the exchange of quark and antiquark in the final state. All structure functions 
(Ji (i = U, L, F) contain the electroweak couplings of the process under consideration. 

For the calculation of the different terms in Eq. (|1]) (corresponding to Born and 
virtual corrections as shown in Fig. la,b) it is more convenient to rewrite the differential 
cross section in terms of all possible parity-parity combinations i,j = V,A that occur 
through the virtual 7 and Z states: 
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i,j=V,A 



(2) 



where 6 and if are the the polar and azimuthal angles of the scattered quark. Here, L 
and H denote the lepton and hadron tensor, respectively, and as usual Na = 3 takes 
into account that the produced quark pair comes in three colors. Because of cylindrical 
symmetry around the beam axis (as indicated in Fig. 2) the integration over (f is simplest. 



The mass parameter v = yl — Am? /q^ (m quark mass, q total energy-momentum transfer) 
enters Eq. (|^) via the two-body phase-space factor PS2 = f /Stt. Note that the product of 
lepton and hadron tensor with fixed parities L^^ Hjj^ is multiplied by the corresponding 
coefficients g^^ 

g^^ = Ql-2QqVeVqRexz + {yl + al)vl\xz?, (3) 
5,^^ = -Qqa^aqReXz + '2veaeVqaq\xz?, (4) 
= {vl + al) al\xA\ (5) 

which contain the neutral-current couplings Vf = 2T/ — AQjsin^ 9^ and af = 2T/, and 
Qf denotes the fractional charge of the fermion. The Breit-Wigner form of the massive 
gauge boson is characterized by 

„ i\/r2 2 

X.{q')= , r ^7. = -^- 4.299- 10-^ GeV-l (6) 

- Mj + iMz Tz 8v2 vra 

The normalization of the symmetric and antisymmetric components in the tensor 
decomposition 

L^" = L^^'^) + Lt^'^] (7) 
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has to be in agreement with Eq. (0) 

L^^^^ = l(p^^p':+p-p^i-igVO, (8) 

r 1 4? 

LiH = (9) 

where p_ and p+ refer to the momenta of the electron and positron, respectively, and 
g = p_ + p+ is the total momentum transfer as shown in Fig. lb. 

In Fig. 2, we show the kinematic configuration of the quark-antiquark final state 
in the center-of-momentum system (cms). The initial e~^e~ beam is aligned along the 
2;-axis, so that the polar angle 6 coincides with the scattering angle between the outgoing 
quark momentum pi and the incoming electron momentum p_. As usual, the square 
of the cms energy is given by = (pi + ^2)^ = {P- + Recall that the remaining 

degree of freedom is captured by the rotation angle ip around the beam axis. In this 
coordinate frame, the particle momenta take a particularly simple contravariant form and 
the orientation of the quark momentum is given hy p = sin 6 cos (p, sin 6 sin cos 6 



Notice the meaning of f = y 1 — Am? /q^ as a dimensionless scale of the quark's velocity 
in the cms frame. 

In the on-shell renormalization scheme the self-energy corrections vanish so that only 
the Feynman diagrams Figs, la and lb are required to calculate the hadron tensor for the 
complete 0(a) process with two quarks q{pi) and q{p2) in the final state. We choose the 
Feynman gauge for the internal gluon propagator. Then, the general substitutions made 
to arrive at the massive QCD vertices from the Born term expressions are simply 

7„ - (l + -4)7„-S^^!i^, (10) 

W ^ (l + C)7„7!, + g '^'^J"''' T5, (11) 

with the induction of charge and magnetic moment form factors in the vector vertex and 
an axial-charge form factor in the axial vertex. 

The explicit real parts of the chromomagnetic form factors are given by 
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m ([^) + 2) ( - m A + ln(l - v^) - 21n2^ 

+F(v 



'12) 
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A-K V \1 + V ) 

ReC = ReA-2Re5, 



(13) 
(14) 



where A is a cut-off parameter to control the soft infrared divergences. At one-loop order, 
the cutoff A is related to the usual dimensional regulator e = A — N hj the following 



correspondence rule [TQ] 



In A 



- 7b + In ' 



(15) 



4-A^ V . 

where 7^ is the Euler constant and ^ the 't Hooft mass associated to the couplings in N 
dimensions. Furthermore, the shorthand F{v) was introduced 



F{v) 



Lis! 



2v 



Lio 



2v 



TC 



+ 3v\n(^^]-4. (16) 



V — 1 J \v + 1^ 

As usual Cp = 4/3 is the Casimir operator defining the adjoint representation of the 
SU(3)c color group. Note that the identity C = A — 2 B connects the gluonic vertex 
corrections of the vector and axial-vector currents for massive fermions. 



The form factors Eqs. ([T^)-(|Tj) have been derived in Ref. [§] with a completely anti- 
commuting Dirac 75 in dimensional regularization following the prescription of Chanowitz 
et al. |TT|. The same final expressions |]7|,P^ are also obtained by employing dimensional 
reduction. In this particular method, the Clifford algebra of the Dirac matrices is reduced 
to 4 dimensions while space-time is still kept in = 4 — e dimensions to regularize the 
otherwise divergent loop integrals. At first glance, naive dimensional reduction and di- 
mensional regularization seem to differ by a finite term. However, at one-loop order the 
inclusion of a global counterterm effectively restores the reduction in the spin degrees of 
freedom, so that both schemes are consistent with each other [|T^. 

In dimensional reduction, a straightforward calculation now gives the 0{a) virtual 
corrections to the Born approximation. The key ingredients of Eq. (0) are 



L^>"''^ Hj^^^ivirtual) = Aq^ [{1 + 2Re A) [2 - v'^{l - cos^O)} 



L^^"^ H^^]""^ {virtual) 
L^^""^ Hf^^^ivirtual) 



+2ReBv^(l-cos^t 



16g^(l + ReA + ReC)w cos^, 
4 g2 (1 + 2 ReC) t;2(l + cos^^). 



(17) 
(18) 
(19) 



The soft IR divergences emerge in the QCD form factors A and C as the cut-off 
parameter A — 0, whereas B gives only a finite correction. The imaginary part ilmD 
in H(^j;) can not contribute to the cross section and finally disappears because of leptonic 
current conservation. Note that the contractions of the symmetric lepton and hadron 
tensors Eqs. (p^ and (|l9|) contribute to the C-even structure functions cr^.^, whereas the 
contraction of the asymmetric tensors Eq. ([T8|) yields ap which is C-odd. 



At this point, we stress that in the explicit 0(a) calculation of the C-odd contribution 
ap, one faces immediately the well-known problem of how to extend the four- dimensional 
75-matrix to dimensions, or, in other words, how to treat the Lorentz indices of the 
totally skew-symmetric tensor e^'^'''^ in A^ dimensions. Neither does dimensional regu- 
larization nor dimensional reduction avoid a direct confrontation with the 75-problem: 
in both cases space-time is A^-dimensional. Moreover, the global reduction counterterm 



of Ref. [|1J] had been taken directly from the massless vector-current vertex. While its 
unrestricted usage for loops with even 75 has been demonstrated, there still remains to 
investigate its significance for odd 75 calculations. 

On the other side, there exists a well-established technique based on the following 



substitution of the axial- vector current 13,14 



7m75 ^ ^5^e/'^^7p7<.7r (20) 

with 

Z5 = l- —Cp + Gia). (21) 

TT 

The finite renormalization Z^, is needed in addition to the conventional wavefunction 



renormalization Z2 to restore anticommutativity [15| which is evidently violated after 
making the replacement Eq. (^). Although this replacement method has been developed 
in first place for massless multiloop calculations (see e.g. ||l6l), the generalization to the 
massive case is straightforward, since the renormalization of the axial current is related 
with the ultraviolet sector, whereas the coUinear and soft divergences relate to the distinct 
infrared sector. 

To complement our calculation, we examine the effects that this particular regular- 
ization scheme has on the gluonic correction of the axial-vector current. For this purpose. 



we consider the contraction of the epsilon tensor with the chiral vacuum polarization ten- 
sor (with cuts in the two final fermion lines of the corresponding diagram). This gives 
the following definition of the C form factor: 
27r 

^ = (g2)2^2 e(/i, i^,Pi,P2) (22) 

d^k Tr (/i + m) 7^ (j/i + ^ + m) ^^^^ (-/a + ^ + m) 7" (/a - m) 7^ 
(2vr)^ A;2{(pi + A;)2-m2}{(p2-A;)2-m2} 

where A; is the running momentum of the gluon loop. 

Then, we can show that within the framework of the 75-replacement prescription 
Eq. (pOD the previous C form factor Eq. (p!^ is recovered 
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2 (1±^ In (1±|) - 2) ( In - 1 ln(l - .^) + In 2 + 1) 



-4t;ln (y-^) +i^(^) +4 



(23) 



In Appendix A, we outline the derivation of C in Eq. (^) with a special emphasis on 
the renormalization procedure and explicitly point out the differences to the dimensional 
reduction approach. 

To obtain the result (^31) we have used the 0(a) expression for the renormalization 
factor Z5. The a^-term in Eq. (|2TD is required to restore the usual chiral Ward identity 
corresponding to the conservation of the 0{a) axial current in the fermionic zero-mass 
limit. 



With the 75-prescription of Eq. (pOD, we obtain the following expression for the 
C-odd part in the differential cross section 

L^^"^ Hl^]^^ {virtual) = 4(4 - 3e) q\l + Re A + Re C)v cos 6, (24) 

where the product of the two epsilon tensors has been replaced by the following natural 
extension of the four- dimensional identity to dimensions 

e(/i, a, /3)e(/i, z/, p, a) -2 5ap V - S^Jisp ) • (25) 

Note that the e-term in Eq. (^) gives additional finite contributions in the virtual part 
due to the infrared poles in A and C when compared with the result Eq. (|1^) obtained 
by dimensional reduction. 
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3 Real-Gluon Emission 



The 0(a) tree-graph contributions of Fig. 1(c) are required to cancel the IR/M divergences 
in the virtual parts, and thus render the physical cross section finite in concordance with 
the Kinoshita-Lee-Nauenberg theorem |]T7|. In the following, we distinguish two types of 
IR/M singularities that occur at one-loop level in the real and virtual parts. 

The soft divergences result from the massless character of gluon field and are reg- 
ulated by introducing a small gluon mass nig. The corresponding regulator A = m?g/q^ 
is compatible with the cut-off definition in the form factors Eqs. (12) and (|lj), and de- 
fines the scale to discriminate the soft-gluon from the hard-gluon domain. In the case of 
real-gluon emission, A 7^ has the specific task to slightly deform the critical boundary 
region of the notorious three-body phase-space. On the other hand, the coUinear diver- 
gences emerge when the gluon field couples with another massless field, and they typically 
manifest themselves as singularities in the mass parameter ^ = Aim? /q^. 

Deriving the hadron tensors /i^^ = V,A) for the 0(a) tree-graph contributions, 
we only have to consider the lowest-order results in the gluon-mass expansion since the 
soft divergences, which will occur after performing the phase-space integration, are at 
most logarithmical, and thus vanish in the limit A ^ when multiplied by an additional 
A. Four-dimensional trace algebra gives the following explicit results for the parity-even 
contributions 



Svr a, Cp 



+4 



f 4 f 4 2 
^ + - + ^ + --2 

y z yz 



-2^ 


-2'- 
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y 
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\z2 


yz 



yz 



(26) 



r^AA 



^TT a, Cf 



g^u < 4^— + ^—ir- + 4^— + ^- 



+4 
+4 



-2 

Pit^Pi. 



y y 
(1-0(2-0 



yz 



2^-2^ 

y. 



z 



Q y y yz 

g2 yz 



g2 1 yz 



(27) 



where the phase-space variables 



y 



2pi- q 



and 



2p2- q 



[28) 



q- 

specify the energies carried by quark and antiquark in the cms frame. Energy-momentum 
conservation in combination with current conservation eliminated the gluon momentum 
P3 and lead to drastic simplifications in these final expressions. 

In the computation of the C-odd contribution these favorable features are absent. 
Especially the 75-replacement prescription Eq. (|^) yields quite lengthy expressions in the 
original trace 



t^VA 



-|Tr(j/i- 



ra) 



1/1 + 1^3 + m , -^2 -I^3 + m 
7a — — — 7[m + 1\m — — 7o 
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-^2 - I^3 + m 



2^2- P3 

7^^175 + 7i.]75 7 



(29) 



2p2-p3 2prP3 
At this point, we shall therefore not reproduce the full expressions for the uncontracted 
hadron tensor but postpone further discussion of the C-odd part until we will come 
to the final contracted results. 

The qqg phase-space is considerably more complicated than the previously discussed 
two-particle complement. In general, five parameters are needed to completely character- 
ize a specific configuration in three-particle phase-space. The phase-space parametrization 
is crucial for obtaining closed analytical solutions after integration out all variables except 
the cos 6'-dependence and the cms energy. Our particular choice for these five parameters 
is ^/q^, y, z, cos 9, and ip. 

Energy-momentum conservation confines the momenta of the three outgoing parti- 
cles to within a plane. Fig. 3 illustrates the angular orientation of the {pi, P2, p^j-pla^ne in 
the cms coordinate frame. The vector normal to the outgoing 3-jet plane is 



n 



— cos 9 cos (fi, — cos 9 sin 93, sin 6' ^ , 



(30) 



and the direction of the scattered antiquark p2 is obtained by a rotation ^ of p = Pi/\pi\ 
around n with the angle x- We find 



P2 = '^nix) P = P cos X — n X p sin X- 



m) 
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Now it is straightforward to write all final momenta in terms of the chosen phase-space 
parameters (with the standard metric g'^'^ = diag(l; —1, —1, —1)) 



Pi 

P2 
P3 



\4^{y + z) (l; Ps), 



(32) 



where the gluon is scattered off in the direction 



P3 = -P^ — P2- 



(33) 



y + z y + z 

Note that p^ = automatically requires to have unit length which yields in combination 
with Eq. (|1|) 

y + z + yz + ^-1 

cos X = P • P2 = r. . . r— (34) 



;i-i/)^-ev/(i-^)'-e 

which makes evident that x gives no additional degree of freedom but is fully determined 
by energy-momentum conservation. It is now easy to obtain the following relations to 
classify the angular dependence that the various hadronic components of Eqs. (|26D and 
(p7|) give upon contraction with the lepton tensor 



L'-'g,. = -4, 



P2tiP2. 



j^^u Plf^P2u + PluP2^, 



^l-cos^Xcos^^^) {(1-^)'-^}, 
\siTi^e [y + z + yz + i - ly 



(35) 
(36) 

(37) 
(38) 



Hence, we readily find the following contracted results after substituting for cos x accord- 
ing to Eq. O) 



vv 



4(1 + 0(2-0- 



yz 



2cos^ 



1 



- 2(1 - 0(2 - 30- - 01 - OS - 4(1 - 0(3 - 20- 



-01 - if- - 201 - 04 -V- + ^z + 2- + 2(7 - 50- 

z'^ z^ z^ y z 
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2 1 

+8^ + 2(2 - 0(1 - 0'— + 2^ + % + 2yz 

z yz z 



(39) 



and 



^TT a, Ce 



8(1-0) (- + -)- 201-0 [h 



4(1 + 



^ + f 



1 2cos6' 

+4 1-0 2-0- + 7^ ^— 

yz [l-yy-^ 



e-2(l-0(2-30--Ol-0'4 

y y 



-4(1 - 0(3 - 20- - 01 - 0'4 + 2^(1 - 04 + 2(1 - 0- - + 4(1 - 0^ 

Z Z Z Z Z 



z 



.y 



^y' , ... .^2 1 



+2(1 - 0- + 2(7 - 60- - 4(2 - 0- + 2(2 - 0(1 - - + 4?/ - A^y' 
y z z yz 

+2{l -i)yz\ 



(40) 



Out of the five variables that parametrize the three-particle phase-space we only need 
to remove the y- and z-dependence by integration to yield the differential cross section, 
depending solely on Ecms = Vq^ and cos 6*. Note that in Eq. (|]) the y9-integration has 
already been included in the two-particle phase-space factor PS2 apart from the flux 
factor. Thus, the appropriate {y, 2;)-integration can be put into the form 



WUreal) 



PS, 



"'fJ.L' 



IGvr^ V 



y+ z+[y) 
dy 



(41) 
for the 



y- z-{y) 

which gives the final real-gluon contributions that have to be included in Eq. 
cancellation of all IR divergences with the corresponding virtual parts. 

The complicated upper and lower bounds of the nested integral in Eq. (^) are 
given in Appendix B where we discuss in more detail the qqg phase-space. In our scheme, 
the phase-space integrals are grouped into distinct classes {Jj}, {§«}, and {di}- Each 
class refers to a different functional dependence on the quark's velocity in the three-body 
system. The individual elements of {Jj}, {§i}, and {di} are identified in Appendix B 
which also treats some specific properties of the integrals and their interconnection. 

Thus, the final integrated results for the 0(a) tree graphs with VV and AA parity- 
parity combinations are 



2^ Ail/ ]^VV/ r\ '^s 



(real) = Cf 

2V TT 



16(1 + 0^2 - 401 + 0:^3 + 8J4 + 4(1 + 0(2 - 0^5 
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+2 cos^^^ - 8^2 - 5^1 - 2(1 - e)(2 - 3e)a2 - e(l - O'as + + 2de + ^Sr 



-4(1 - 0(1 - 2oa8 - 01 - 0% + 2(2 - 0(1 - O'aio + 2aii - 2(1 + 5oai2 



+201 - oai3 - ai4 + 2ai5 



(42) 



se^i - 16(1 - e)a2 - 401 - 0J3 + 8(1 + 0:14 



+4(1 - 0(2 - 0^5 + 2 cos^^ I - 4(2 - 0^2 - eai - 2(1 - 0(2 - 30a2 - 01 - 0% 

+4a4 - 4ea5 + 2(1 - oae + 4(1 - oa? - 4(1 - o'as - 01 - o'ag + 2(2 - o(i - o'a 
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+2(1 - oaii - 2(1 + 2oai2 + 201 - oai3 - au + 2(1 - oa 



15 



(43) 



Note that {Ji}, {§«}, and {di} are the smallest units of the three-particle phase-space 
one arrives at after parametrization in the energy variables {y, z) and subsequent partial 
fractioning. These single components are process- independent. With the additional in- 
tegral class 7i (see Ref. 0) they constitute a complete set of solutions for any one- loop 
real-emission process from massive fermions. Differential or total observables such as pro- 
duction rates or various polarization components can all be expressed succinctly in terms 
of these units once they are derived. 

To contract the asymmetric part of the lepton tensor L\^^'^^ with the 75-odd hadronic 
complement Eq. (|29|) , we use the following identities 



= 1 2v cos fc', 



e(/i, z/,pi,p2) 



L 



—i 2v cos 6, 

2cos^ 



(44) 
(45) 



{(2 -102/ +(1-10^ 



j^^^u <f^, Pu P2, q) P2u - ejjy, pi,p2, q) P2^ ^ ivcosOy, 



(^2)2 

e(/^,Pi,P2,g) Qu - e{iy,Pi,P2,q)q^ 



2^2 



iv cos 9 {y + z). 



(46) 
(47) 

(48) 



In these expressions, the contraction of two epsilon tensors with one common index has 
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been replaced by an appropriately antisymmetrized product of metric tensors which is 
well-defined in either four or dimensions: 



e n ,rMPo"^ — cAp n'^^ 



(49) 



Then, the scalar product of two Lorentz vectors is invariant in different space-time di- 
mensions so that Eqs. give unique results. 

It is now straightforward to obtain the following contracted results for the general 
case of A^- dimensional Clifford algebra 



= 2,271 as Cp cose (l - |e) 

-2(4-30j + ef + 252/ + (2 + 5)^+{2-.(l-ie)}^ 
+ {6-5(l-|e)}--(2-^)^ + (l-0(2-0- + 0(^ 



y W ^ J 



(50) 



Putting e in Eq. ( pH ) to zero and then performing the (y, 2;)-integration gives the 
following result for the C-odd tree-graph contribution within dimensional reduction 



L^^ Eltireal) = L^^ E^^ {real) 

■(4-50S2-^(l-0(S3 + §5) 



o 9^ a 
2 Lp cosO 

V n 



2(4 - 30§4 



+^§6 + 2§8 + 2§9 + 6Sio - 2§n + 2(1 - 0(2 - OS 
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(51) 



On the other hand, adopting the 75-prescription of Eq. ( PP] ) we obtain 



L>^-^H^;^{real) = L^'^H^^{real) 



V TT 



e^Si - (4 - 50§2 - 01 - 0(§3 + S5) - 2(4 - 3OS4 



+^§6 + (2 + e)§s + {2 - £ (1 - } §9 + {6 - £ (1 - |e) } § 



-(2 - £)Sn + 2(1 - 0(2 - 0§i2 + 2£Si3 + 0{e] 
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(52) 
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2 Cp cosO 

V n 



(4 - 50§2 - (1 - - 0(S3 + §5) - 2(4 - 30§4 



+^S6 + 2S8 + 2S9 + 6S10 - 2Sn + 2(1 - |£)(1 - 0(2 - 0§i2 + 0{e) . (53) 

Eq. (l5^ ) displays the full result neglecting 0{e) bracket terms, which do not contribute 
in the limit e ^ 0. Since the logarithmic divergences within the massive gluon scheme 
correspond at 0(a) to the 1/e poles of a dimensionally regularized theory, we can not 
naively discard all 0(e) -terms. In Eq. (|53|) only the non- vanishing 0(e) bracket terms are 
retained. However, these additional 0(e) terms represent the finite difference in the 75- 
odd real-gluon correction between dimensional reduction and dimensional regularization 
with the 75-replacement scheme. Finally, when the corresponding virtual-gluon correction 
Eq. (p^ ) is added, these contributions disappear, and one obtains the same total result as 
in dimensional reduction Eq. (^Tj). 



4 Complete 0{cx) Differential Cross Sections 

With the explicit expressions for the virtual and real contributions of the previous sections, 
we are now in the position to give the full analytical results for the 0(a) differential cross 
section of the annihilation process e~^e~' —^'~f,Z—>-qq. Full quark- mass dependence was 
kept and no other approximations were made so that these analytical formulae are valid in 
a closed form over the entire energy range of perturbative QCD. All angular dependence 
resulted naturally from basic kinematics of the two- and three-body final state. 

A severe test on the total results is provided by the cancellation of the soft and 
coUinear divergences separately for each angular distribution within the corresponding 
real and virtual parity-parity combinations. Using the relative phase-space factor Eq. ( PT| ) 
in the sum of virtual and real gluon parts, we obtain the following total 0(a) results 
(including the Born level) 



L^'^H^J' (total) =4q' 



;i + 2 Re i) {2 - t;2(l - cos^^)} + 2 Re 5 v^{l - cos^ 



all 

+ 7^ - - 8(1 + 0^2 - 2^(1 + 0^3 + 4^4 + 2(1 + 0(2 - 0:^5 



An V 
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+ cos'^^ { - - 5eai - 2(1 - 0(2 - 3oa2 - e(i - ifh + + 2a6 + 4a7 
-4(1 - 0(1 - 2oa8 - 01 - O'ag + 2(2 - 0(1 - O'aio + 2aii - 2(1 + 5oai2 



+20i-oai3-eai4 + 2ai5}| 



(54) 



Lf"" H^,"^ (total) = 



;i + 2ReC') t;^(l + cos^ 



a 1 \ 

- 4^^! - 8(1 - o:J2 - 201 - 0^3 + 4(1 + 0^4 + 2(1 - 0(2 - 0:^5 

47r V [ 

+ cos^e { - 4(2 - 0^2 - eai - 2(1 - 0(2 - 3oa2 - 01 - O'as + 4a4 - 4^05 

+2(1 - Oae + 4(1 - oar - 4(1 - O'as - 01 - O'ag + 2(2 - 0(1 " 0% 



+2(1 - oaii - 2(1 + 2oai2 + 201 - oai3 - + 2(1 - oais} 



(55) 



where the wiggle on top of form factors or integrals denotes quantities free from soft 
divergences as A — 0, which essentially means that only linear terms in In A have been 
discarded. We find that in the massless limit ^ — * all coUinear singularities cancel, as 
expected. It is now easy to obtain the full 0(a) expressions for the structure functions in 
Eq. (|l]) by using the convolutions 



cr, 



a \ 



ij=VV,AA 



(56) 



where the unpolarized and longitudinal projectors are explicitly given by 

+1 



Jdcos9 {5cos2^- l}, 

-1 
+1 

JdcosO {2-5cos20}. 



(57) 
(58) 



Similarly, one finds that the remaining total VA combination is IR finite. Adding 
Eqs. (p^ ) and (|53D gives for the C-odd component of the differential rate 



a 



'l + ReA + ReC)v 



a 1 \ ~ ~ 

+7^ - - (4 - 50S2 - 01 - 0(S3 + §5) - 2(4 - 30§4 + ^§6 
47r V 
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+2(S8 + §9 + 3Sio - Sii) + 2(1 - 0(2 - OS12} , (59) 

where we substituted A and C according to Eqs. ([T^ ) and (p3[) without soft divergences. 
Note that for the A^-dimensional 75-replacement scheme the terms proportional to e in 
the real part have exactly canceled with the extra contributions in the virtual part. Thus, 
dimensional reduction and dimensional regularization give the same result for the forward- 
backward asymmetry in the differential cross section. The independent result for dimen- 
sional reduction is easily obtained by adding Eqs. (|^) and (^T]). 

In the massless calculation, there are no QCD one-loop contributions to the forward- 
backward asymmetry p|,p2|. Apart from the cancellation of the spurious IR divergences 
in Eq. (p9D, we recover this specific result in the limit ^ 0. On the other hand, we have 
checked explictly that integrating the VV and AA parts of the differential rate over cos 6 
reproduces exactly the analytical results for the total cross section as found in Ref. [Q. 

The numerical estimates of the (D{a) differential cross section for bottom quark 
production with a fixed-point mass mi,{mh) = 4.3 GeV are shown in Fig. 4. The running 
of the strong coupling is implemented by taking af'^M^) = 0.123 in the modified minimal 
subtraction scheme (MS) for five active flavors. For the bottom quark at energies above the 
next flavor threshold and in general for the top quark, we apply the appropriate matching 
conditions for six active flavors using the corresponding one-loop QCD renormalization 



group equations [1 



Fig. 4a gives a surface plot of the 0(a) differential rate as a function of the cms 
energy Ecms = \f(? and the cosine of the scattering angle cos^^. On the Z-peak, the 
asymmetry is clearly pronounced and yields for E^ms = = 91.178 GeV values ranging 
from 5800 pb {6 = tt) to 8300 pb {6 = 0). The minimum is located at cos 6 ~ 0.16 with 
3494 pb. For higher energies the differential rate falls rapidly off to give at 100 GeV cross 
sections of the order of only 100 pb. Nevertheless, the massive 0(a) corrections become 
increasingly important off the Z-peak. In Fig. 4b, the 0(a) results are compared with the 
Born approximation. Note that the dominant 0(a) contributions are below the Z-peak 
in the domain < cos^ < 0.4 and above the Z-peak in the domain —0.5 < cos 6 < 0. 
Thus, the area of maximum correction (4.3-4.5 %) consists of a strip of approximately 7r/6 
width that shifts from the upper to the lower hemisphere when passing the Z-threshold 
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(see also Fig. 2). At 100 GeV the corrections amount for cos 6* = 0.3 to full 4.5%. 

Fig. 5a depicts a similar three-dimensional plot with the differential cross sec- 
tion for e+e" — > 7, Z — > at one-loop QCD level. The value for the top mass is 
mt{mt) = 174 GeV. Note that we have chosen an energy range sufficiently higher than 
the tt threshold to avoid full interference with the non-perturbative sector. The domi- 
nant contribution to the differential rate is given by top quarks scattered in the forward 
direction. Along the forward direction, we find at ~ 425 GeV a peak value of ~ 0.85 pb 
for the cross section which includes an 0{a) correction of nearly 20%. At 350 GeV the 
strong corrections have an impact of 52 %. Thus, it is conceivable that for top quark 
production non-perturbative effects still prevail in energy domains considerably above the 
tt threshold. 

In Fig. 5b, we give the distribution of the energy and scattering regions most signifi- 
cant to the strong corrections. Note that the 0{a) terms contribute dominantly along the 
forward-scattering axis where the differential rate takes maximum values. On the other 
hand, bottom production yielded the most important 0{a) corrections perpendicular to 
the beam axis which corresponds to minimal values on the saddle surface in Fig. 4a. 

5 Conclusions 

In the present work, we have derived from first principles the full analjd;ical 0(a) results 
for the differential cross sections in heavy-quark production. No mass approximations or 
further restrictive assumptions were made. All angular distributions arose naturally from 
the underlying phase-space kinematics. 

The systematic treatment of the tree-graph contributions relies on the exact integral 
solutions of the massive three-body phase-space. These integrals are divided into several 
classes according to their functional dependence on the final quark velocity. A complete 
collection of these phase-space integrals allows to describe any one-loop bremsstrahlung 
process including mass effects. 

In particular, we put emphasis on the consistent treatment of the axial- vector cur- 
rent to handle the spurious chiral anomalies that are usually present in dimensionally 
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regularized calculations with an odd number of 75 's. We explicitly showed that dimen- 
sional reduction and dimensional regularization with the 75-replacement prescription yield 
identical results for the forward-backward asymmetry. 

In our numerical estimates, we found that at one-loop level the final state QCD 
corrections already modify the Born approximation by approximately 4.5 % (100 GeV) 
for bottom production and more than 25 % (375 GeV) for top production. We presented 
details on the differential distribution, and made the essential observation that for the 
bottom quark the 0{a) corrections are dominant in the scattering region perpendicular 
to the beam axis whereas for top quarks 0{a) contributions become very important along 
the forward-scattering axis. 

We conclude that full analytical calculations for heavy-fermion pair production at 
one-loop order are feasible and provide an attractive alternative to existing Monte-Carlo 
generators. 

Acknowledgements. M.M.T. wishes to thank Sergei Larin for clarifying comments on 
the treatment of the axial-vector current within dimensional regularization, and Sergio 
Novaes for a stimulating discussion on the numerical results. He further gratefully ac- 
knowledges the support given by the Alexander-von-Humboldt Foundation. This work 
has been supported by CICYT under grant AEN93-0234 and IVEl under grant 03-007. 
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Appendix A: Renormalization of the Axial- Vector 

Current 



In this appendix, we give a brief account of the subtleties associated with the renormal- 
ization of the axial-vector current at QCD one-loop level. 

As we explained in Section 2, dimensional regularization with a conventional anti- 
commuting 75 or a replaced 75 produces identical results for the gluonic vertex correction 
to the axial- vector current. Explicitly, we gave meaning to the 75-matrix in a dimension 
other than four by the replacement rule 

7^75 Z5-e/'"^7p7^7^, 

where the finite renormalization constant is necessary to reinstate the validity of the 
axial-vector Ward identities in the final results [|T3|,n^ . 



Using Eq. (^Qf ) one can recast Eq. (p2D into the following expression for the unrenor- 
malized form factor 

d^k Tr (/i + m) 7^ + ^ + m) 7p7^7^ {-^2 + ^ + m) 7° (/a - m) 



X 



(2vr)^ ;t2{(pi + A;)2-m2}{(p2-A;)2-m2} 

where all indices are A^-dimensional. Note that the normalization has to be in agreement 
with the overall A^-dependent factor in the Born contribution Eq. (|2^). 

In the solution of this one-loop integral we use the following structure for the two- 
and three-point functions [pf = p\ = m^) 

h{Pl,P2) - 7 (2^)^p|(p^ + ^)2_^2}{(p^_^)2_^2} - ^"^'^ 



I^{Pl^P2] 



d^k ki" 



(27r)^ p {(pi + ky - m2} {{p2 - ky - m^} 
Q(v)fil^. (A3) 



d^ k k^^ k^ 

^'^^P'^P'^ = I {27r)^ k^ {(p, + ky - m^} {{p, - ky - m^} ^^^^ 
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1 



^M^^r(pi>P2)- (A5) 



Notice that the last identity only holds if one takes g^'^ g^u = N = 4 — e. All coefficient 
functions P, . . . ,T are dimensionless and depend only on the mass parameter v. 

An explicit calculation yields for Eq. ( |Al| ) 

1 + 



1 - H + 2Q(^)} - (l + i^) Riv) 



2v 



+ (l - H Y^^^M - 1(4 - e)T{v) + 0{e] 



(A6) 



Here, the e-dependence stems exclusively from the A^-dimensional trace algebra so that 
the result for dimensional reduction is obtained from the above formula by putting e = 0. 
We verified this result also in an independent calculation. 



Adopting the conventions of Refs. [13,20|, we can further decompose the elements 



of Eqs. 



into the so-called Passarino-Veltman functions 



P{v) 
Q{v) 
R{v) 
S{v) 
T(v) 





m 


ijj, ^ 1 




L 


(47r)2 








l-i;2 " 


4(47r)2 






i/i"^ 1 




L 


(47r)2 16t;2 




i/i ^ 


Bo (m' 


4(47r)2 



12 2 2 

[q ; m , m 



m 



^2. 2 2 



m 



Bo (g^;m^m2) - 2 



(A7) 
(A8) 
(A9) 
(AlO) 
(All) 



For conciseness, we do not consider the imaginary contributions to the form factor but 
concentrate on the real parts of the following general scalar integral representations 

1 



Bo [pl;ml,ml^ = {i 



2\-l / jN 



rk 



|A;2 _ 77^2 II (^/j; + pi)2 — mfj 



(A12) 



[it: 



2\-l / jN 



rk 



|/c2 _ j72^||(/;; _)_ p^)2 _ ^2| _|_ _ ^2| 



(A13) 
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In the required Passarino-Veltman functions we neglect terms proportional to e and 
drop the imaginary parts. These scalar two- and three-point functions are divergent as 
e ^ and read 



Bq (m^; 0, m 



ID ( 2 2 2 



- - 7, + In ( ^ ) + 2, 
2 



Co {rn^ 1 q^i m?] 0, m^, m^^ 



1 



- 7^ + In f ) + 2 + In {^—^ 

\ J V 1 + f 

3t;2 1 



(A14) 
(A15) 



V 



£ V / 1 + f ^ 



In 



1 -t; 
l + v 



(A16) 



where -F(f) was previously defined in Eq. (|T6|). Here, we have included the 't Hooft mass 
/i which naturally arises when the couplings are extended to dimensions. 

Thus, the additional real contribution induced by e in Eq. is given by 



5C '. Ch„n C 



'IGni ttgCp e R{v) 



a, 
2^' 



(A17) 



Note that the first term in brackets on the right-hand side depends on the mass parameter 
V and is thus connected with the infrared structure of the theory. On the other hand, the 
remaining term stems from the ultraviolet sector. 



The UV divergences in Eq. ( |Al| ) are removed by renormalization. Apart from the 
conventional on-shell renormalization scheme we take into account the renormalization of 
the axial- vector current particular to the special 75-prescription used [|T3|. Therefore, the 
entire renormalization program amounts to isolating C^ep in the following equality 



1 + Crcp — Z2 (^1 + Cba 



(A18) 



where the subscript refers to the particular 75-definition chosen. In the MS scheme, the 
quark-field renormalization constant Z2 (with gluon mass regularization) and the finite 
axial-vector renormalization constant at one-loop order are given by 



Z2 



a 



2 , fiTlfi^' 

h 7i3 - In ' 



4-21nA + 21n^-41n2 



(A19) 
(A20) 
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where the e in generates spurious finite contributions when multiphed with the IR pole 
terms in the virtual- and real-gluon parts. However, these finite terms cancel in the total 
results. 

In the intermediate step of the derivation, we get after eliminating the UV diver- 
gences with Z2 the relation 

^2(l + C,..,) = l + C + ^C, + 5C, (A21) 

where C is the chromomagnetic form factor for dimensional regularization with anticom- 
muting 75, namely Eq. (p^. By multiplying the right-hand side of Eq. ( [A21| ) with Z5, we 
finally obtain the fully renormalized result (neglecting imaginary parts) 

Crop = C + 5C — —Cp 
In 



2 1^1±^ In {j^^ - 2^ In A5 - i ln(l - + In 2 + 1 j 



-4i;ln(^^j +F{v) + A 



which agrees with the result Eqs. (O) and (p3| 



(A22) 
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Appendix B: Properties of qqg Phase-Space Integrals 



The massive three-particle phase-space is as usual defined by 



PS, 



3 f]N 



(2n 



,N+3 



(Bl) 



where in this particular case of fermion pair production mi = m2 = m. Further, we put 
^3 = and = 4. Regularizing by a small gluon mass has the advantage to shift 

all complications to the boundary functions of the integrals whereas dimensional regular- 
ization produces an e dependent integration measure which does in general not permit to 
fully exploit powerful substitution techniques. It appears that the most complicated of 
these phase-space integrals are only solvable in four dimensions. 

Integrating over the gluon momentum and all the internal angles except for x gives 

1 



PS, 



(An) 



d\P^\ \Pi\^ 

1=1 



J dx sinx 5((g -pi - P2) 



Ipf + m? 

where the delta function expresses energy-momentum conservation 

4 



(B2) 



5 -P2)' 



X 



cosx 



yz + y + z + i-l + 2k 



V 



(B3) 

Note that we do not put A = in the argument of the delta function which directly effects 
the deformation of the phase-space boundaries to regulate the IR singularities. 

Finally, we obtain with the (y, 2;)-parametrization given by Eq. ( pH]) for the three- 
body phase-space 



PS, 



y+ z+(y) 



dy / dz 



(B4) 



y- z-{y) 

where the upper and lower bounds of the nested integral follow directly from the constraint 
-1 < cosx < +1 and Eq. (B3) 

y~ 
z±{y) 



(B5) 



21/ 
4y + e 



A 1 



2 y yV V 
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In the calculation of the differential rates one requires the solutions of the following 
three distinct types of phase-space integrals 

dydz fi{y,z), §i = [ ^ fi{y,z), 



, i-y)'-^ 

'''' (B6) 



where fi{y,z) are real rational functions in the quark-energy variables y = 1 — pi - q/ 
and z = 1 — p2 - q/ q^ ■ Apart from the numerical verification of each individual analytical 
solution there exist several consistency checks among the different integral classes. For 
brevity, we will only discuss relations among {Jj} and {Sj}, but similar considerations 
apply to 

It is important to recognize that although {§j} and {^j} emerge as new integrals (in 
addition to {Jj}) in the transition from the the total cross section to the differential form 
of the cross section, this procedure can not fundamentally alter the divergence structure 
in the soft divergences Consider the following argument: The integrand f{y,z) of a 
divergent § is first integrated over z to obtain F{y). Then, before executing the remaining 
y-integration, we add and substract a suitable function F{y): 



dy 



jdy m+J^ 



F{y)-J{l-yy-iF{y) 



(B7) 



Now F{y) can be choosen that way that the second integral on the right-hand side of 
Eq. (p7|) becomes regular as A ^ and all soft divergences are contained in the first 
much simpler integral of type J. It becomes clear that § and J exhibit the same divergence 
structure in A, i.e. for any ^ 7^ 

S(ln) = J(ln) -|- finite terms as . (B8) 

On the other hand, the coUinear singularities provide another useful tool to establish a 
correspondence between S- and 3-integrals. For A 7^ 0, we can in general find a suitable 
combination of §j and with i ^ k so that 

Si(ln{, 1/0 - Si(ln^, 1/0 = Ji(lnO 1/0 + finite terms as { ^ 0. (B9) 
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Furthermore, simple algebraic manipulations yield valuable identities that save consid- 
erable labor, e.g. the following substitution was useful in the derivation of the real-tree 
graph contributions 

An exhaustive list with all three-body phase-space integrals relevant to the deriva- 
tion of the 0{a) differential cross section for massive fermion pair production is given 
below. The integral sets {3i} and {Sj} with the exception of §13 have been pubhshed be- 
fore 1^. The integral class {3i} is entirely new. Some of its components are also important 
in the angular dependence of the quark's alignment polarization [IS]. For the ^-integrals 



we use the shorthand w = y (1 — v^) /(I + v^). The standard book on dilogarithms and 
their integral representations is Ref. |pT| . 

Class Integrals 



"Ji = J dydz 



\v (1 + |e) - |e (1 - \i) In 



[ dydz — = ( dydz — 
J y J z 



1 + v 



J dydz — = J dydz — 

-y (inA^ + In^ -21nw -21n2 + 1) + 2 - In (^^ 



dydz - = dydz - 
J z J y 

.1, (5 - i{) + 1 (1 + le) In (ii^) 
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^5 — dydz 



yz 



= (-21nA5 -ln^ + 41n?; + 21n2)ln(^^ 



+2 



Li2 ( ^— r-) - Li2 ^"'^ 



+ 3 



2w 
1 -i; 



-I 



Class S Integrals 



Si 



o?y dz 



V(i-y)^-^ 



So — 



dydz 



= 2 In 



-lnA5 + lln^ + ln(l- y^) -21n(2- y^) +ln2- 1 



dz 



2; 



^ + LI2 



2 Li, 



2 y ■ - - V 2 - ' ' 

ln(l + ln(l - v) 



dy dz 



27 



In - 1 In e + ln(l - V^) - In (^) + In 2 



r dydz y_ 
4 



/ 



dydz 

z 



1 



32 



i2-(2 + er-^±^e' + 2(8-oe i^^zTi 



-|lne + ln (2-^ 



' -2 



z 



+ (2-lf)to(H-^)-V^+2„ln 



l-v 
1 + v 



28 



>11 



/ 



dy dz 



^(1-1/)^ 



= y- 



Li2 ( -4—) + Li2 '^^ 



-ln(i±^)ln 



l-v 



+ 3f In 



-2Li.(iVe) + iln^(K) 
) + 1(18 + 0- 1(20 -OV? 



l-v 



1 + t; 



>12 



dydz 



^ In ( 2 In + I In ^ + 4 ln(2 - y^) - 4 In ^; - 4 In 2 - 2 In (^^ 



1 - V 



V V\/e(2-x/e). 



^2 /ve(Wl)V /Mliiiv^S). 



2 

+ - 

V 



l-v 
1 + t; 



+ Lis 



(1 



,V?(2-Ve). 



+- 



T- /1+^^ 

Ll2 



Li2 



1 - tj 



2?; 



1 -t; 



Lin 



2v 
1 + v 



>13 



dy dz 



1 

16 



-e In e + 2 ln(2 - ^) + 4(2 - - 4 (2 - C^) 



Class Integrals 



ai = 



(i-i/)^-e 
i-e, /i + t' 



. In , 
4-e Vl-^; 



r dydz 1 
y {l-yy-i y 
6 /l+^^\ 



29 



/ 

4 



dy dz 1 

InA^ -ln^ + 21nv + 21n2- 1 



24 , /1 + v 
m 



e(4-0 Vl-^; 



/ 



dy dz 



(i-y)'-e 



6 \, + 
m 



4-a \i-v 



y 



2 + k + ie-^)ln^'^'^ 



A-i \l-v 



-1(6-0^ 



/ 



dy dz 



yz 



(i-y)'-^ 



12 



e (4 - 0^ 



In 



1 + v 



/ 



dz 



(i-y)'-e 



3e /i_^Vnri±^u 



4-n vi-^;; 4-e 



as = 



dy dz 1 



1 

7! 



Li2(w) - U2{-w) + Li2 K ^ ^ ?w 



a. = / 



dz 1 

{i-yy-^ ^ 



30 



4 



lnA5 -ln^ + 21n't; + 21n2-^-r-^ln'^'^"^^ 



2v 



1-v 



r dydz 1 



(r 


+ v\ 




- vj 



-41nA2 -|lne + 51n(l + V?)+41n(l- Vf) 



1-e 



Li2 ( ^^^r-) - Li2 '^^ 



ln(2 + V{) + 61n2 
3 



+ 



1-e 



Li2 



2?; 

1 -V 



2v 
1 + v 



+ 



1 



1 



Li2(w) — Li2{-w) + Ll2 I 7= w I - Ll2 I 7^ w 



1-e 



Lio 



l + w 



Li2(^)+Li2((2 + V^) 



1 + 



-Li2((2 + Ve)^^KLi2 



2v^ 



,(2 + Ve)(l + «^). 



- Li2 



2v^ 



,(2 + Ve)(l-«^), 



/ 



dy dz 



(i-y)'-^ y 

12 24 

-1 + 



4-e (4-0^ 



In 



l + v 



2v 
4^ 



/ 



dz y 



(i-y)'-e ^ 

i InC + ln(2 + - ln(l + v^) - 2 ln2 



^ VI — u 



1- 



Li2(w) - Ll2(-W) + Ll2 K _ ^ w I - Ll2 I _ ^ w 



+Li2 ( ^) - Li2 { + Li2 f (2 + ^ ) - Li2 ( (2 + Ve) 



1 — w 



+Li2 



2Ve 



,(2 + v/e)(l + «^). 



2Ve 



,(2 + Ve)(l-^), 
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)13 



/ 



dy dz y 



dl4 — 



dy dz 



Z^ 



In 



1 + v 
1-v 



2v 



a 



r dy dz y^ 

J (l-y)2_. - 
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i(3+e)in 



l+v 
1-v 



i In e + ln(2 + - ln(l + V^) - 2 In 2 



(1 - V^Y 



Li2(w) - Li2{-w) + Li2 [ ^ ^ w 



Lis 



2 + 

2-Ve 



+ (3 + 



1 + w 



+ Li2 (2 + x/e)— — 



1 — -u;' 

-Li2((2 + +Li2 



40-16£ + £2 n + v 
H r+ — —In 



2v^ 
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1 -i; 



V(2 + ^/e)(l + w)^ 
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Li2 



2v^ 



,(2 + Ve)(l-^) 
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Figure Captions 



Fig. 1: Feynman diagrams contributing to da/dcosO {e~^e~ — > qq) up to order (D{a): 
(a) Born term, (b) virtual corrections, and (c) gluon bremsstrahlung graphs. 
Diagrams (b) and (c) give the definition of the particle momenta. 

Fig. 2: Two-particle kinematics for the process e"'"(p_|_) e~(p_) — > q{pi) q{p2) in the 
cms coordinate frame. The e+e~ beam line coincides with the z-axis. The 
scattering angle between the electron momentum p_ and the quark momentum 
Pi is 9, and angle (f defines the orientation around the beam axis. 

Fig. 3: Three-particle kinematics for the process e"'"(j9+) e^(p_) ^ q{pi) q{p2) diPs) in 
the cms frame. The momentum definitions agree with those of Fig. 2 except 
for the additional gluon momentum p^. The vectors (pi, p2, Ps) span the 
qqg production plane. The angle x = ^{Pi,Pi) a-nd the quark momentum pi 
uniquely specify all the remaining momenta within the plane. 

Fig. 4: Differential cross section for bottom quark production: (a) 0{a) result as 
a function of the cms energy and cos^, and (b) the 0{a) corrections com- 
pared to the Born approximation. Shown are the values for the ratio 
1^ da{0{a))/ dcos 6 / da{BoTn) / dcos 6 — 1 in percent. The bottom quark mass 
is mb{mb) = 4.3 GeV. 

Fig. 5: 0{a) Differential cross section for top quark production with a top mass of 
nT'ti'iTT't) — 174 GeV: (a) surface plot to show functional dependence on the cms 
energy and scattering angle, and (b) the 0(a) corrections compared to the 
Born approximation. At 375 GeV the corrections amount to 30 %. 
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